Ordered phases of XXZ-symmetric spin-1/2 zigzag ladder 
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Using bosonization approach, we derive an effective low-energy theory for XXZ-symmetric spin- 
1/2 zigzag ladders and discuss its phase diagram by a variational approach. A spin nematic phase 
emerges in a wide part of the phase diagram, either critical or massive. Possible crossovers between 
the spontaneously dimerized and spin nematic phases are discussed, and the topological excitations 
in all phases identified. 
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I. INTRODUCTION 



The appearance of unconventional spin-liquid phases 
in frustrated Heisenberg models around a critical point 
separating different magnetic ordered phases is a long 
standing, intriguing issue which has attracted notable 
theoretical and experimental interest in recent yearsi. 
From a theoretical point of view, this is quite a chal- 
lenging problem which calls for a deep reexamination of 
standard theories. Conventionally, the stability of a spin- 
ordered phase may be investigated by spin-wave theory or 
by more sophisticated field-theoretical approaches based 
on the non-linear a-model. Spin-wave theory is in prin- 
ciple able to detect instabilities at any wave-vector, al- 
though a reliable description would require a systematic 
1/S expansion, S being the magnitude of the spin. The 
non-linear tr-model offers a better description of the crit- 
ical behavior close to an instability point, yet it has a 
major limitation. Namely, it takes into account only the 
long-wavelength Goldstone modes of the ordered state 
under consideration, but fails to describe the excitations 
at the wave- vector of the competing ordered state. Hence 
it does not provide any information about the phase 
which can emerge under increasing the effect of frustra- 
tion. 

One-dimensional spin models might turn useful to at- 
tempt an improvement of the field-theoretical approach 
in view of an extension to higher dimensions. The sim- 
plest example of one-dimensional frustrated spin model 
is the spin-1/2 Heisenberg chain with antifcrromagnetic 
nearest neighbor exchange J ' and frustrating next near- 
est neighbor exchange J (for a recent review and refer- 
ences therein, see Ref 2). Besides the general interest, 
this model is also relevant for realistic materials, such as 
CsiCuCl^, in which magnetic Cu ions get arranged into 
a zig-zag fashioni 3 . 

In the classical limit, the J-J' spin chain has a Neel 
long range order for j = J/J' < 1/4, with characteristic 
momentum q = tt/clq. In the Neel phase time reversal 
symmetry is preserved if combined with a translation by 
one lattice spacing ao. For larger values of j, a spiral 
ordering at momentum satisfying cos (g do) = — 1/4 j is 



stabilized at the classical level. There parity and time 
reversal symmetries are separately broken. 

Quantum fluctuations modify the classical phase 
diagram 4 . By Mermin- Wagner's theorem, spin rota- 
tional symmetry cannot be broken in one dimension: 
the Neel long-range ordered phase turns into a quasi- 
long range ordered one characterized by power-law de- 
caying correlations of the staggered magnetization. The 
low-energy effective critical theory is the level- 1 Wess- 
Zumino-Novikov-Witten (WZNW) model (free massless 
bosons with central charge C = 1). On the contrary, the 
spiral order disappears completely in favor of a sponta- 
neoulsy dimerized phase. The transition point is sligthly 
shifted with respect to the classical value, j c ~ 0.241^. 
Within the WZNW model formalism, the transition is 
driven by a perturbation which is marginally irrelevant 
(relevant) at j < j c (J > j c ). At j = j c a Bcrczinskii- 
Kosterlitz-Thouless transition takes place, and an ex- 
ponentially small spectral gap opens up in the region 
j > j c 7 . The system continuously passes to a two-fold de- 
generate, spontaneously dimerized, massive phase. Upon 
further increasing j, the gap reaches its maximum at the 
exactly solvable Majumdar-Ghosh points, j = 1/2, after 
which it slowly decreases^. Even though the ground state 
of the J-J' chain remains dimerized at all j > j c , above 
the Majumdar-Ghosh point the system reveals signatures 
of the classical spiral phase: the spin-spin correlations 
beco me incommensurate, as it was shown numerically in 
Refs l7l9t Since this occurs far away from the region of 
applicability of the SU(2)i WZNW model, there is little 
scope to improve the field-theoretical description of the 
gapless phase at j < j c to account for incommensurate 
correlations that emerge well above j c . More promis- 
ing is to approach this problem from the opposite side, 
j ^> 1. If J' = 0, the even and odd sublattices of the 
spin chain decouple, and the model effectively describes 
two decoupled Heisenberg chains, one for each sublattice. 
Classically this corresponds to the case when the spiral 
wave number is equal to q — 7r/2ao. Switching on a small 
J ' transforms the model to a weakly coupled two-chain 
zigzag spin ladder, with the interchain coupling giving 
rise to a marginally relevant perturbation that opens up 
a gap and brings the system back to the dimerized phase. 
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In addition it should also move the relevant momentum 
q away from 7r/2ao towards 7r/a . Therefore incommen- 
suration and the spectral gap are supposed to appear 
together in this limit, which makes a field-theoretical de- 
scription more plausible. 

Indeed, in the limit j 3> 1, a novel, par ity-breaking 
(twist) perturbation was identified in Ref.[l(Jas a natural 
sourse of the spin incommensurabilities. The twist term 
has a tendency to support a finite spin current along the 
chains which would account for the expected shift of the 
momentum. However, for the SU(2)-symmetric zigzag 
ladder, the situation still remains rather unclear. Ap- 
parently, the appearance of a nonzero spin current is not 
compatible with the requirement of unbroken spin rota- 
tional symmetry (see, however, the discussion in sec. VII). 
On the other hand, no reliable information about the ac- 
tual role of the twist operator at the strong-coupling fixed 
point can be extracted from the Renormalization Group 
(RG) analysis^ because of the perturbative nature of this 
approach. Thus, the structure of the low-energy effective 
field theory for the SU(2)-symmetric S=l/2 zigzag ladder 
still remains unknown. 

The situation changes much to the better in the pres- 
ence of strong spin anisotropy. Namely, close to the XX 
limit, a self-consistent, symmetry-preserving mean-field 
approach shows that the twist operator can stabilize a 
new, spin-nematic (chiral) phased. In this doubly de- 
generate phase a nonzero spin current polarized along 
the easy axis flows along the ladder, and the transverse 
spin-spin correlations are incommensurate and may even 
decay algebraically within some parameter range, as it 
is the case at the XX point. This picture is supported 
by recent numerical simulationsiii2ii2iI4. In particular, 
the numerical work by Hikihara et aiii*i£ has indeed con- 
firmed the existence of the critical spin-nematic phase in 
a broad region of the phase diagram for spin-anisotropic 
chains, both for integer and half-integer spins. NMR 
experiments with compound CaV^O^^ having a spin- 
1 double-chain zigzag magnetic structure have revealed 
the gapless nature of the spectrum, which may be an in- 
dication to the critical chiral state of this material. On 
the other hand, numerical simulations have revealed the 
existence of a new gapped chiral phase in a very narrow 
region between the dimerized phase and the critical chi- 
ral phase, but, within numerical accuracy, only for inte- 
ger spins. In that gapped phase, the spin current coexists 
with dimerization; accordingly, the spin-spin correlations 
are incommensurate but decay exponentially. Recently 
the phase diagram for general spin S has been studied 
by bosonization technique^ and by means of the non- 
linear cr-modelii, verifying the existence of both critical 
and gapped chiral phased for integer spin. 

In this paper we present a detailed study of the 
phase diagram of the XXZ-symmctric, frustrated, spin- 
1/2 chain in the limit J' <C J. Using a variational anal- 
ysis of the bosonized Hamiltonian we identify possible 
phases of the model. In addition to the critical spin- 
nematic phase and to the commensurate spontaneously 



dimerized one, we find conditions for the existence of a 
massive spin-nematic region for the S — 1/2 case. We 
also characterize the topological excitations which occur 
in each region of the phase diagram. 

In the following section we introduce the model and 
discuss the bosonization approach. In section III we 
demonstrate how the variational approach can be applied 
to the twistless ladder, in which only the dimerization op- 
erator plays a role. Critical spin nematic phase, driven 
only by the twist operator, is studied in section IV. The 
interplay between dimerization and twist operator and 
other emerging phases is discussed in section V. In sec- 
tion VI we discuss a ferromagnetic phase which turns out 
to be dual to the critical spin nematic phase. In section 
VII the RG approach is implemented to study the in- 
terplay between different twist operators at the border 
of these mutually dual phases. The last section contains 
conclusions. 



II. THE MODEL AND ITS LOW-ENERGY 
LIMIT 

We consider a frustrated spin- 1/2 Heisenbcrg chain 
with 2L sites, described by the Hamiltonian 

2L 

H = [Ja &n &n+l + S£ S^ +2 ] , (1) 

a=x,y,z n—1 

where 

J x = J y = J > 0, J z = JA, 
J' x = J y = J'> 0, J' z = J' A'. (2) 

In what follows, A and A' will be treated as independent 
anisotropy parameters. Upon the transformation 

S2n — > Si(n), S 2n+ i — > S 2 (n), 

the model JTJ is mapped onto the zig-zag spin-1/2 ladder 
Hamiltonian 

L 

H = E EE J a S?{n)S?{n+l) 

a=x,y,z n—1 i= 1,2 

+ E E J ^i>) + Si>+i)]s 2 >) (3) 

a—x^y,z n 

Let us discuss some general properties of this model. 
In the limit J = 0, describes a standard Heisenberg 
antiferromagnetic chain where the spin-spin correlations 
are modulated with wavevector q = tt. On the contrary, 
when J' = 0, the equivalent model J3J describes two 
decoupled spin chains. The modulating wavevector in 
this case is q — ir/2. When both J and J' are finite, we 
may expect two possible behaviors of the spin structure 
factor S(q) = {2L)- 1 E„,m< S « ' S m )e^"-™) : either it 
is peaked at q = ir and q = tt/2, or it shows a single 
peak at an incommensurate qo which smoothly moves 
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from q = tt/2 at J » J' towards q = tt when J' 3> J. 
Translated into the zig-zag ladder language, the former 
case implies that, for n large, 

(Si(n) • Si(0)) = (S 2 (n) • S 2 (0)) 

= F (|2n|) + F n (\2n\) + (-1)" F T/2 (\n\), 
(Sx(n) • S 2 (0)) = F^\2n + 1|) - #(|2n + 1|), 
(S 2 (n) • S x (0)) = F '(|2n - 1|) - F^\2n - 1|), 

where F (|n|), F ff (|n|) and F„./ 2 (|n|), as well as the 
primed ones, are smooth real functions describing the 
contributions of the q = 0, q = n and q = tt/2 modes, 
respectively. The difference between the 1-2 and 2-1 spin- 
spin correlators, as well as the absence of inversion sym- 
metry n — > —n, reflect the fact that the zigzag ladder 
lacks two Z 2 symmetries - the 1 2 interchange symme- 
try and site parity Ps (understood as Pg ' <giPg). How- 
ever, if the model is gapless or possesses a small spectral 
gap inducing a macroscopically large correlation length, 
then site-parity is effectively restored at long distances. 

If, apart from q = 0, the spin structure factor has a 
peak at an incommensurate wave vector qo £ [7r/2,7r], 
then we expect that 

(Sx(n) • Si(0)) = <S 2 (n) • S 2 (0)) 

= F (\2n\) + F qo (\2n\) cos2ng , 

(Sx(n)-S 2 (n))=i^(|2n+l|) 

+ ^ o (|2n + l|) cosg (2n + l), (4) 

(S 2 (0)-S 1 (n)) = ^(|2n-l|) 

+ ^ o (|2n-l|)cosg (2n-l). (5) 

We notice that the presence of the modulating factors in 
®,0 makes the the breakdown of Ps even more pro- 
nounced and, contrary to the commensurate case, this 
breakdown will survive the continuum limit we are going 
to adopt. Thus, the two different types of spin correla- 
tions - commensurate or incommensurate - can be dis- 
tinguished within a continuum, low-energy description by 
an asymptotic restoration or breakdown of the site-parity 
symmetry. 

As discussed in the Introduction, in this paper we are 
going to study the model or equivalently JJ), in the 
limit J 3> J' of weakly coupled chains. That allows us 
to adopt the well-known continuum description of each 
XXZ chain based on the bosonization approach^ (see 
also Ref. for a recent review) and then treat the in- 
terchain coupling as a weak perturbation. Bosonization 
of the XXZ zigzag spin- 1/2 ladder has already been dis- 
cussed in Ref|2JJ. Here we review this procedure in more 
detail paying attention to the structure of the effective 
continuum model which is important for the subsequent 
analysis of the phase diagram. 

We start with the Abelian bosonization of a single XXZ 
spin-1/2 chain. Its universal, low-energy properties in 
the gapless Luttinger liquid phase (— 1 < A < 1) are ad- 
equately described by the Gaussian model for a massless 



scalar field ip(x) = <pr{x) + <Pl{x), 



Hxxz —> J dx Hg(x), 

H G = (v s /2) Iq- 1 (d x tp) 2 + Q {d x df 



(6) 



Here = —<pr(x) + <Pl(x) is the field dual to y>(x), 
n(x) = —d x "d(x) being the momentum conjugate to 
(p(x), v s ~ J a,Q is the velocity of the collective spin 
excitations, and Q is the the Luttinger liquid parame- 
ter which determines the compactification radius of the 
field ip. The dependence Q = Q(A) in the whole range 
— 1 < A < 1 is known from the Bethe-ansatz solution: 
1/Q = 1 — (1/tt) arccos A. Thus Q varies in the range 
oo > Q > 1 when A takes values within the interval 
— 1 < A < 1. In particular, Q = 2 at the XX point 
(A = 0), and Q — 1 at the SU(2)-symmetric (Heisenberg) 
point (A = 1). Throughout this paper it will be assumed 
that A < 1 (i.e. Q > 1). In this case the perturbation 
to the Gaussian model JSJ, Xu cos \/87r(p (Xu ~ JA), 
which in terms of the Jordan- Wigner fermions originates 
from Umklapp processes, is strongly irrelevant and will 
be dropped in what follows. 

Since for each chain only the uniform and staggered 
low-energy modes survive the continuum limit, the cor- 
responding spin densities can be parametrized as follows: 

Si(ra) — > a Si(x), (x = na ) 

Si(x) = 3 itR (x)+3 i!L (x) + (-irni{x) = 1,2). 

Here ao is the lattice spacing, 3i.R,L are chiral compo- 
nents of the smooth part of the magnetization of the i- 
th chain, and is the staggered magnetization. The 
latter is even under site parity transformation (Ps) and 
odd under link parity transformation (Pl)- A distinc- 
tive feature of the zigzag ladder is that it is invariant 
under mixed parity: P^ ®P^ and P^ ®P S 2 \ By this 
symmetry, strongly relevant terms, n^n^, which deter- 
mine the spin-liquid properties of the unfrustrated spin- 
1/2 ladders^, are instead forbidden in model ©■ As a 
result, in the low-energy limit, the interchain perturba- 
tion, H' = Jijj + Wtwist ■ is contributed by the "current- 
current" interaction^^ 



n 



J J 



a—x,y : z 



n T a T a 
ya<Jl <J2 ; 



and also by "twist" termsifi allowed by the P, 
symmetry: 



(i) 



n 



twist 



= \ ^2 g a at)T a + ^-g a T , 



(7) 



>(2) 



(8) 



Here 



T a = n\d x nl-nld x nl, T° = e x d x e 2 - e 2 d x e u (9) 

are chirally asymmetric operators with conformal spin 1, 
and ti ~ (— l) n Sj(n) ■ S,(n + 1) represent the continuum 
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limit of the dimerization operators. (Note that for a sin- 
gle chain e(x) is even under P^ and odd under Pg.) The 
coupling constants are given by g x — g y = g± = J'ciq, 
9z = g\\ = J'A'a Q . 

The "current-current" and twist perturbations are of 
different nature. The former are parity [i.e. Pg^Pg^] 
symmetric. If acting alone, provided that the interchain 
exchange is antiferromagnetic, these lead to spontaneous 
dimerization of the ground state (see section III), the 
existence of massive topological excitations (spinons), 
and the onset of short-ranged commensurate interchain 
spin correlations^*^. The twist terms, whose appearance 
stems from the frustrated nature of interchain interac- 
tion, explicitly break parity. However, by the previous 
discussion, parity can be broken either in a mild way, 
which is the case when the leading asymptotics of spin- 
spin correlations are still commensurate, or more pro- 
foundly, i.e. explicitly inducing incommensurations in 
the spin correlations. Both patterns of the low-energy 
behavior of the system will be discussed below. 

In model only the vector part of the twist pertur- 
bation, g a T a , emerges in the continuum limit. The scalar 
part, goT°, although absent in the bare Hamiltonian, is 
generated in the course of RG flowiS. 

For this reason we will assume that such a term is 
present at the outset, with a bare amplitude go . 

Let us first bosonize Hjj. In terms of the rescaled 
fields, fa = (l/y/Q)ifi and 9i = y/Q'&i, the "currents" 
Ji=J?,R + JZ L ,™e given bj^ 



C 



R = R + = l/R_ which considerably simplifies the per- 
turbative analysis. 

Performing an additional rescaling of the fields, fat = 
\/R±$±, 0± — y/RzfQ±, for the transverse (xx, yy) part 
of Hjj one finds: 



Hjj., ± = 2g ± V J?J* = -± (V + T) , (13) 
* — ' 7ra 

a—x^y 

V = cos ^AttK + <& + cos y/4-!rK-Q-, 
T = cos cos 



where Aj^ = g±( 2 /ira, and 

K+ = QR, K_ = R/Q. 



(14) 



(15) 



To bosonize the twist perturbation we use the 
bosonization formulas for the staggered magnetization of 
the S=l/2 XXZ chain (see e.g. Gj|): 



nl = - (C z /na) sin \J1-kQ fa 
nf = (CJira) exp(±iy/2n/Q Ot), 



(16) 



where C a (Q) (a = x, z) are noniniversal parameters 
(thei r exac t dependence on Q was recently found in 
Refs. l25l26|) . Then, in terms of the fields $±, 9±, the 
twist term becomes: 



H 



LwisL 



= E Xi(Di ' 



(17) 



i=l,2,3 



e l V QZ cosy/2nQfa, (10) with 



where a is the short-distance cutoff of the bosonic the- 
ory, and C(Q) is a nonuniversal (and yet unknown) pos- 
itive constant approaching the value 1 in the SU(2) 
limit. Using the definitions (|10fl and passing to the 
symmetric and antisymmetric combinations of the fields, 
fat = {fa ± fa)/V2, 9± = (6>i ± 6 2 )/V2, we find that the 
longitudinal (zz) part of Hjj adds to the Gaussian part 
of the model transforming the latter into 

H g - E T WW + K\d^f] , (11) 



er = ± 



with 



1 

r± 



— = A 1± 



.911 Q 



irv s 



1 ± 



9\\Q 
2irv s 



0(gh 



(12) 



The exact dependence of R± on the dimensionless pa- 
rameter g»Q /ttv s is unknown. Therefore we will restrict 
ourselves to the case \9\\\Q/ttv s <C 1 and keep only lin- 
ear terms in the expansion (|12fl . For a weak interchain 
interaction (|<?|||/7n; s -C 1), this is justified almost for 
the whole range |A| < 1 except for a narrow region 
A + 1 ~ (g\t/wv s ) 2 close to the ferromagnetic transi- 
tion point. The parameters R± then satisfy the relation 



A i ~ C x9±/ a , A 2 , 3 ~ C 2 z {±g\\ + go) /a, 



(18) 



and three bosonized twist operators 0i,2,3 related to 
T* (i = 0,1,2,3) as follows: 




d x S+ sin y/AirK-e-, (19) 



-in 



(20) 



d x <$>- sin y/4jrK+<f>+. (21) 



Thus, the bosonized continuum version of our model, 
H = Ho + Hjj-± + Htwistj represents a Gaussian field 
theory of two scalar fields, 



Hn 



E n » ] = E t fa-*-> a + ^ 



(22) 



with perturbations and (fT7|) which couple the (+) 
and (— ) channels together. Since Hq is perturbed in a 
relevant way, the relationship between the coupling con- 
stants of the original model JTJ and the parameters of H , 
obtained within our weak-coupling approach, is not to be 
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trusted. For this reason we will consider those parame- 
ters as independent. Namely, we will treat the Hamilto- 
nian H, as a low-energy effective theory for a most gen- 
eral class of frustrated zigzag spin- 1/2 ladders, sharing 
the same symmetry properties with the model Q). 

The scaling dimensions of the perturbing operators are: 
d v = K + + K_, dr = K_ + -^, 

di = l + lf_, d 2 = l + -^, d 3 = l + K+.{23) 

K— 

Their relevance (d < 2) or irrelevance (d > 2) can be un- 
derstood from Fig.l where the plane (K-, K + ) is shown. 
The point K + = K- = 1 corresponds to the SU(2)- 
symmetric zig-zag ladder where all perturbations (includ- 
ing the Umklapp term) are marginalia. This point and 
its close vicinity will not be discussed in this paper. Due 
to the condition K + /K_ — Q 2 > 1, the physical part of 
the (if-, K + ) plane lies above the line K + = K_ and can 
be divided into four sectors in which at least one twist 
operator is relevant: 

sector A : d\ < 2, dp > 2, d<z, d 3 , d-p > 2, 

sector B : di < dp < 2, cfe, d%, d? > 2, 

sector C : dp < c?i < d 3 < 2, d%,dj: > 2, 

sector D : d 2 < 2, d T > 2, di, d 3 , d-p > 2. (24) 



1 





\® 

® X Ox 


® 

o 2 




Ox V 






© 




Ox Oi v/ 





FIG. 1: The parameter space of the model 

Notice that, except for the operator T>, all other per- 
turbing operators have a nonzero Lorentz spin: Si, 2, 3 = 
1, = 2. Strictly speaking, the conventional criterium 
of relevance does not apply to such operators (see e.g. 
Ref. Il9t) because in higher orders of perturbation the- 
ory they can generate relevant scalar perturbations. Us- 
ing standard fusion rules for the Gaussian model, we 



have analyzed the structure of various terms appear- 
ing in the second order of perturbation theory. There 
are marginal terms leading to small corrections to the 
parameters K± and the velocities v±, as well as those 
which rcnormalize the already existing coupling con- 
stants. Besides, new scalar operators cos 2y / ATrK + ^> + , 
cos2 v /4tt/X_$_, cos 3^4nK-0- and cos 2^4^X16- 
are generated. The first two of them have scaling di- 
mensions 4QR and AQ/R, respectively, and are therefore 
strongly irrelevant (since Q > 1, R ~ 1). The third 
operator has dimension 9K^ — 9R/Q and becomes rel- 
evant roughly at Q > 9/2, which is a region of the fer- 
romagnetic intrachain exchange, far away from the XX 
point. This region will not be considered here. Finally, 
the last perturbation ~ cos 2-y/47ri^_6_ has dimension 
4_ftT_ = 4R/Q and thus becomes relevant at A"_ < 1/2, 
which corresponds to a vicinity of the XX point Q = 2. 
However, even in that case its role is subdominant, as we 
will show later. 

Thus, the continuum model we will be dealing with in 
the remainder of this paper reads: 

H = H Q - (XJira)V+ E \ t O % . (25) 

i=l,2,3 

where Ho is given by (|22ll . For later purposes, we have 
suitably inverted the sign of the coupling constant Aj^ by 
making a shift of the field <!>+: $ + — > $ + + ^tt/AK + . 
In what follows, we will analyze possible phases of this 
model in the four sectors A,B,C,D by using a general- 
ization of the standard variational approach^! that ac- 
counts for ground states with nonzero values of topolog- 
ical charges, (d x Q-) and (d x $>±). The very possibility 
to incorporate such states within the variational method 
stems from the fact that the twist operators (|19|) - H21|) are 
products of fields belonging to different Gaussian models 



III. TWISTLESS LADDER 

Before addressing the role of the twist terms in 1251) it 
is instructive first to apply the variational approach to 
a simpler frustrated two-leg ladder model 23 - 28 which, in 
the continuum limit, is free from parity-breaking pertur- 
bations yet being spontaneously dimerized. This model 
is the two leg-ladder version of the standard J1-J2 frus- 
trated Heisenberg plane, in which the interchain coupling 
includes besides the usual on-rung coupling, J±, a frus- 
trating exchange, J x , across the diagonals of the plaque- 
ttes. In the XXZ case its Hamiltonian reads: 

a.n ?'—1.2 a, n 

+ E J'a [S?(n)S$(n + 1) + Sf(n + l)S 3 »(n)] ,(26) 
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where J a are defined as in J3J), and 

J' z = J ± A', 
J" = JxA'. 



J x j y 



The 

lne S(L) 



T" — T" — J 

u x y " X ? 



s ^ reflection symmetry of the model H26|) 
forbids the marginal twist perturbations to appear in the 
continuum limit. The additional condition J± = 2J X 
eliminates the nln^ part of the interchain coupling and 
thus makes the two decoupled Gaussian models only per- 
turbed by the current-current interchain interaction Q 
with coupling constants g± = 2J x A'ao and g\\ = 2J x ao- 
So the bosonized continuum Hamiltonian has the struc- 
ture of Eq.J2El with A_i_ 7^ and Ai,2,3 = 0. As we shall 
see below, the dimerized phase it describes also occurs in 
the most part of sector C of the zigzag ladder, Ea. (|2*51l . 
where T> is the most relevant perturbation to the Gaus- 
sian models lj2*2*)l . 

To implement Coleman's variational procedure^!, we 
introduce a trial ground state, 



vac 



) = |0;to + ,to_) = |0;to + ) (g) |0;m_), 



(27) 



which describes free bosons in the (±) sectors with 
masses m±. These are regarded as variational param- 
eters. To estimate the variational ground-state energy 
density, E^m+^mJ) — (vac|7i|vac), one needs to nor- 
mal order the Hamiltonian with the prescription that, in 
the normal- mode expansions of $>±{x) and Q±(x) 1 m± 
should be treated as infrared regulator masses. 
Upon normal ordering^ 



iLq — -'v TO± /Tq 
&ttv± ( a ) 



O 



(28) 



cos y/4irK + $+ = 
cos = 



\m + \a 



K- 



cos y/iirK+Q-i 
cos 



where Af m± are the normal ordering symbols. Subtract- 
ing from l|28|l the diverging contribution of the zero-point 
motion when a — > 0, ignoring for simplicity the difference 
between the velocities v± and defining the dimensionless 
quantities, 

4ira 2 m±a 4Ai a , , 

£= Eo, M± = -±-, z ± = ^^, (29) 



we find that 



S(M+,M-) 



Ml + Ml 



-z±\M 4 



\M- 



(30) 



Withous loss of generality we choose M± to be positive. 
With the condition zj_ <C 1 in mind, it should be under- 
stood that the masses M±, obtained upon minimization 
of £, should satisfy M± <C 1. 



At dp > 2 only a trivial solution exists, M± = 0, cor- 
responding to a critical regime in which the interchain in- 
teraction is irrelevant and the two chains asymptotically 



decouple in the low-energy limit. At d-p = K + 
we find a nontrivial solution, 




k: 



'K+ ^JK 
with ground state energy given by 



2(2-d-p) 2{2-d T) ) 



1-K+ 


)m\ 


2K+ 


2-d D ^ 


)Ml 


2K+ , 



1 - K- 
2K_ 



Mi 



K_<2 
(31) 

(32) 
(33) 



This solution describes a strong-coupling, massive phase 
in which the fields <&+ and 0_ are locked in one of in- 
finitely degenerate minima of the potential U (<£>+, 0-) = 
— (X^/TrajV. Since Aj^ > 0, these minima decouple into 
"even" and "odd" sets: 



-2n A 



4K_ 



-2n^ 



4K, 



(2n+ + 1) , 0_ = 



4K_ 



(2n_ + 1) ,(34) 



where n± = 0, ±1, ±2, .... The existence of these two in- 
equivalent sets reflects two-fold degeneracy of the sponta- 
neously dimerized ground state. Transverse dimerization 
is the order parameter; it is defined as (e±(x)} wher o 7 i 23 



e±(x) = ni(x) ■ ri2(x) oc C 2 cos y/inK-Q 
1 



cos v/47r/X_$_ . (35) 



-~-Cl (cos v^A+$ H 

Since the field 0_ is locked, its dual <&_ is disordered 
and, hence, the expectation value of the last term in l|35(l 
vanishes. Hence, 



(e ± )=±e , e C*\M_\ K - + ^-C'i\M + 



1 



\K 4 



(36) 



with the two signs of e corresponding to the even and 
odd vacua, respectively. 

The discrete (Z-i) symmetry that is spontaneously bro- 
ken in the ground state is generated by even-odd interset 
transitions of the fields, 



A$+ = ±^71-/4^+, A6_ = ± v / ir/4:K-. / 



(37) 



and is related to translations by one lattice spacing on 
one chain only. This is not an exact symmetry of the mi- 
croscopic Hamiltonian (|26|l but rather appears as an im- 
portant property of the corresponding low-energy model 
with a "current-current" perturbation. The excitation 
spectrum of the model consists of pairs of massive topo- 
logical kinks (spinons) interpolating between two adja- 
cent minima of the potential hi. The kinks carry two 
topological quantum numbers - the total spin 



k r 

-± / da; d x $+{x), 



(38) 
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and the relative longitudinal spin current 



dx d x O-(x), (39) 



which, according to (|37|l . take fractional values ±l/2 29 . 



value of the current operator l|41|) which is to be deter- 
mined self-consistently. The variational procedure is the 
same as in the previous section with the exception that 
the ground state energy in the (+) channel will acquire 
a piece proportional to (j+) 2 . Otherwise this sector re- 
mains gapless: M + = 0. Using dimensionless notations, 



IV. CRITICAL SPIN NEMATIC PHASE 

Now we are coming back to the continuum model 125(1 
for the XXZ zig-zag ladder. We begin our discussion with 
sector A where the twist operator 0\ is the only rele- 
vant perturbation to the two decoupled Gaussian models 
H^K Making a shift 9_ -> 6_ + (l/tyy/n/K- we write 
the low-energy model in sector A as follows: 



Ha 




d x Q+ cos y/AwK-B-. (40) 



This model has the same structure as that for the XX 
zigzag ladder considered in Ref[l3 Not surprisingly, the 
variational procedure we will follow now leads to qualita- 
tively the same results as those obtained for the XX case 
within a symmetry-preserving mean-field approach^. 

The interaction term in l|4U|) couples the vertex oper- 
ator in the (-) channel to the topological current den- 
sity d x Q+ in the (+) channel. The latter determines the 
z-component of the spin current which flows along the 
chain direction, j^j, 



3+ = 3\\ = 



-ox/^±d x Q + . 



(41) 



We observe that a finite Ai , see Eq. (14011 , generates an ad- 
ditional contribution to the spin current, which flows 
along the interchain bonds. By the continuity equation 
related to the conservation of the z-component of the 
total spin, one finds that 



3± = 



: Al cos ^jAllK-Q- 



(42) 



The total spin current is therefore j" = j?, and the 

twist operator 0\ is nothing but a coupling term j?, 

The structure of the perturbation in the model 14U|) 
suggests that the ground state admits finite values of the 
mass gap in the (-) channel and the spin current in the 
(+) channel. So one needs to treat both of these two 
quantities as variational parameters. To this end, we 
keep boundary conditions periodic for the field 0_ (x) but 
impose twisted boundary conditions for the field + (x): 



G + (x) = el(x)- 1 -J]ff + x. 



(43) 



Here 0l(x) is a massless harmonic Bose field satisfy- 
ing periodic boundary conditions, and ji_ is the average 



J+ = 



2na 



JK+v"^' V K + v 

for the variational energy density £ we obtain: 



£(J+,M-) = ~ {M 2 _ + Jl) T z±J+M 



(44) 



(45) 



As before, we have chosen Af_ to be positive. The 
(^F) signs in the interaction term correspond to two sets 
of vacuum expectation values of the field 0_: 0_ = 
y/ir/K-n and 0_ = ^n/K^(n + 1/2), respectively. 

Minimizing £ with respect to M± and J7+ we find that 
the (-) channel is gapped, 



M_ 



(46) 



if K- < 1. This is actually the condition d\ < 2 for the 
twist operator 0\ to be a relevant perturbation, which is 
satisfied in sector A. At the same time, the gap supports 
a finite value of the spin current in the (+) channel: 



J+ 



±- 



±K_ z-. 



(47) 



We notice that the dimensionless transverse current 
defined by 



Jl 



2na 



-31 



zi cos 



s/K + u 

also acquires a finite ground-state expectation value 
J± = (J±) =TziM K -, 



(48) 



(49) 



which exactly cancels — J+, so that the total spin cur- 
rent is zero. This results in a spin nematic (or a staggered 
spin-flux) phase characterized by local spin currents cir- 
culating around elementary plaquettes in an alternating 
way. This type of ordering does not break time reversal 
symmetry. In sector A the spin nematic phase is critical 
because the spin-density fluctuations in the (+) channel 
remain gapless. 

We notice that the transverse current can be associated 
with the chirality order parameter. The latter is defined 
as 



(k z (x)), k z (x) = [ni(ir) x n. 2 (x)] 



(50) 



and, according to bosonization rules i|16|) . transforms in 
the continuum limit to 



k z (x) oc cos y/4irK-®^(x) oc j±(x). 



(51) 
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As the dimerized phase discissed in sec. Ill, the spin 
nematic phase is doubly degenerate because the mixed 
parity symmetry P^mPym is spontaneously broken in 
the ground state. The two degenerate phases differ in the 
signs of the longitudinal and transverse currents. Conse- 
quently, apart from the massless bosonic mode describ- 
ing low-energy fluctuations of the total magnetization, 
there exist massive topological Z2 kinks corresponding 
to vacuum- vacuum transitions, 



J+ 



-J+, 0_ -> 0_ ± y/n/AK-, 



and thus carrying the relative spin current jf_/u = ±1/2; 
see Eq. figfy . 

The presence of a finite longitudinal spin current in 
the ground state makes the transverse (xy) spin correla- 
tions incommensurate. Since the (+) channel is massless 
and described by a Gaussian field with a K + dependent 
compactification radius, the correlations will decay alge- 
braically with a nonuniversal exponent. Making use of 
bosonization rules l|16l) . Eq. H43[) an d the fact that the 
field 0_ is locked, one easily finds the asymptotic be- 
haviour of the transverse spin correlation function: 



(_l)x/a 



I1/2K+ 



(52) 



where the wave vector go — tt^7\\/uK + . At the XX point 
(Q = 2, R = 1, K + = 2) the spin correlations decay 
according to the power law |x| 1//2 , in agreement with Ref. 

M 

The ground state energy of the critical spin nematic 
(CSN) phase is given by 



CSN 



1 - K- 
2K- 



M 



(53) 



Before closing this section, we would like to briefly 
discuss the role of the scalar operator cos 2 \/AirK-<d_, 
which is generated in a higher orders and becomes rele- 
vant at K_ < 1/2. In the presence of a finite spin current 
in (+) channel, this term transforms the effective Hamil- 
tonian in (-) channel to a double-frequency sine-Gordon 
model: Hq — A e fj sin yj AnK- _ — g cos 2 ^ AitK_Q_, 
where A e s oc Ai(j'i). It is known^L^i that the g-term 
can induce an Ising transition to a new massive phase if 



l/(2-K_) 
off 



g > and 31/2(1-2*--) > \ 

ity, however, is not satisfied since the amplitude g 
is rather small and, hence, the presence of the second 
harmonics does not qualitatively affect the above results. 



The last inequal- 

A cff 



second twist perturbation, O3, is either irrelevant (as 
in sector B) or the least relevant (as in sector C). In 
Appendix A we explicitly show that its role is indeed 
subdominant in sectors B and C far from the SU(2)- 
symmctric point, K + = K- = 1, 

Thus, the effective Hamiltonian reads: 



H 



B/C 



(+) 



x 



- — cos JAitK + <& + cos WIttKZQ- 
net 

+ -2h = d x S+ sin y/4nK-6- (54) 



The potential in l|54|) contains both the sine and cosine 
of the field 0_; so its vacuum value 01 is expected to 
be located somewhere within the interval (0, ^J-k/AK-) 
and must be such that in a massive phase with M_ 7^ 



(vac I cos v 7 AttK-Q*_ |vac) = M_ 
(vac I sin ^/^kKZQ*_ |vac) = 0. 



(55) 



Setting 0_ = Q*_ — 7/ -J AttK- , we arrive at the following 
expression of the dimensionless variational energy: 



E(J h M + ,M_, 1 ) = ~(Ml 



M 2 _ + Jl) 



T ztJ+M*- sin 7 - z±M* + M*~ cos 7. (56) 

Its minimization with respect to M±, J+ and the angle 
7 yields the following set of equations: 



M K (z±M+ + SU17 =F z\J + COS7) 



0, 



rK 



J+T^iMl sin 7 = 0, 



K 4 



M+(l - z ± K + M+ 



M-(l - z ± K-M* + M_ 



K- 



K- 



COS7) = 0, 

2 

COS 7 
17) = 0. 



(57) 
(58) 
(59) 

(60) 



There are two obvious solutions of these equations in 
which only one of the two perturbing operators is effec- 
tive. In these solutions the angle 7 takes two values: 
and it/2. The corresponding phases are, respectively: (i) 
a fully gapped D phase already described in sec. Ill, with 
zero current (J7+ = 0) and nonzero masses M± given by 
Ea. (|31() . and (ii) a CSN phase with nonzero J+ and M_ 
given by Eqs. gZJ and (gSJ. 

Eqs. H57fl - H6U[) admit one more solution where the com- 
bined effect of the two relevant perturbations leads to an 
intermediate value of the mixing angle 7, 



V. MASSIVE SPIN NEMATIC AND 
DIMERIZED PHASES 

Let us now move to sectors B and C where the prop- 
erties of the systems are determined by the interplay be- 
tween two most relevant perturbations, Oi and T>. The 




cos 7 



and a finite mass gap in the (+) channel, 



M 4 



(61) 



(62) 
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This is a noncritical or massive spin nematic (MSN) 
phase characterized by the coexistence of a reduced spin 
current J + 



J+ = [^+] C SN sin 7 

and a nonzero dimerization 



oc ± 



1 



ClM~- cos 7 + -C 2 Z M + 



(63) 



(64) 



An important observation is that the minimal value 
of the variational energy 1)56(1 is still given by expression 
Therefore, the energies of the MSN and CSN phases 
are related as 



£msn = £csn 



K, - 1 



2K 4 



-Mi 



So in sector B (K + > 1) the MSN phase is energetically 
less favorable than the CSN phase, and the ground state 
should be chosen between CSN and D phases. Accord- 
ingly, in sector C (K + < 1) the competing phases are 
MSN and D. 

Consider first sector B. Here we need to compare the 
ground state energies of the D and CSN phases given by 
Eqs.J2SJ an d 1(53(1 . These are of the same order when 
the mass gaps of the two phases, Eqs. l|^T|) and (g|jj), 
become comparable. Notice that the coupling constants 
Zi and zj_ are both proportional to g±_ and, hence are of 
the same order of magnitude; their ratio is 



where C is a nonuniversal number, 
dition 



(65) 



Therefore, the con- 



z 2 ± - dD ~ zt~ dl (66) 

can be satisfied only in some vicinity of the line K+ = 
1 where scaling dimensions of the operators D and 0\ 
become equal. 

As already mentioned, it is not possible to establish a 
precise relationship between the parameters of the orig- 
inal, microscopic model and the effective low-energy 
theory ((25|) . As a result, the parameter C in 1(651) is un- 
known. Therefore we are forced to consider two cases, 
Z\ > z± and Z\ < z±, on equal footing and draw plausi- 
ble scenarios for each of them, leaving the final choice to 
future numerical work. 

Setting K + = 1 + 5 with \5\ < 1, we find that the 
condition ((66|) translates to the relation 



5 = (1 - 



ln(l/zi) ' 



(67) 



the relation 1(6 7() determines a phase boundary between 
the CSN and D states. The transition is of first order, 
associated with discontinuities of the spin current and 
dimerization order parameters. It can be easily shown 
that in the case Z\ < z±_ the D phase, occupying a narrow 
region close to the line K+ = 1, extends over the whole 
C phase. 

In the opposite case, z\ > z±, Ea. H67|) has no solu- 
tion for 6 > 0, implying that CSN is a stable ground 
state in the whole sector B. Moving to sector C opens a 
possibility for the MSN phase. If 1 — K + <C 1, the condi- 
tion z\ > z±_ admits a small nonzero mass M + given by 
((62(1 . Thus is sector C the upper boundary for the MSN 
phase is K + = 1. The lower boundary is found from the 
requirement cos 7 < 1 (see Eq. I(61() h Within the log- 
arithmic accuracy, this brings us again to Eg. ((67(1 . this 
time for S(K^) < 0. 

Thus, if the ratio Zi/z± > 1, then the CSN and D 
phases are "sandwiched" by the MSN phase occupying a 
narrow region in sector C 



1 -(!-#_) 



hi(zi/z_i_) 
ln(l/zi) 



< K+ < 1 



(68) 



attached to the line K + (see Fig. 2). In all this region 
£msn < £d- The transitions that occur on the upper 
and lower boundaries of the MSN phase are continuous. 
When moving from sector B to sector C through the 
MSN phase, the mixing angle 7 varies from ir/2 to 0. 
Correspondingly, the current J7+ decreases from its nom- 
inal value [J7+]csn and vanishes at the lower boundary, 
whereas the transverse dimerization e±_ increases from 
zero at the upper boundary and reaches its value eo, 
Eg. (1361) . in the pure D phase at the lower boundary (see 
Fig. 2). A possible way of driving the ladder to pass 
through these phases is shown in Fig. 3. 



\CSN 




CSN 






D 

1 / 



K- 



1<- 




D 



"T" 



1<- 



FIG. 2: Possible phase transitions in the ladder: the left one 
is when z± < z±, MSN phase is very narrow, and the right 
one is for z± > z\ 



VI. FERROMAGNETIC PHASE 



This relation determines a line 5 = S(K_) which lies 
entirely in sector B (S > Q) and is located very close to 
the line K + — 1 only if z\ < z±. Under this condition this sector corresponds to the case of a ferromagnetic 



Let us now consider sector D where K + > K- > 1. 
This condition implies that K + K- = R 2 > 1, and so 
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interchain interaction < 0). The effective iow-energy 
model 



Hd=H[ +) +H ( ) +A 2V / ^5 a; $ + sinV47r/A^$_ (69) 

contains only one relevant twist operator and is dual to 
model (I4U|I describing the SN phase in sector A: mapping 
between these two models is achieved by duality trans- 
formations 



2Ai^A 2 , K±->l/K±, $ ± ^6± 



(70) 



Using this correspondence, we can readily translate the 
results of sec. IV to the present case. In particular, the 
spontaneously generated spin current J7+ of the SN phase 
transforms to the z-component of the uniform spin den- 
sity. So the ground state of the system in sector D is 
ferromagnetic (F). Contrary to the spin nematic phase, 
the F phase breaks time reversal invariance but preserves 
parity P| X) ® P ; (2) 



Shifting the field <&_ by */ irK- / 4 and passing to di- 
mensionless notations for the coupling constant, 



Z2 



V 



and total magnetization 



5+ = 



2na 



-m 



we write the variational energy density as 

£ = \{Ml+Sl)Tz 2 S + M 1 J K -. (71) 
Its minimization yields a finite gap in the (-) channel, 



M_=K_ 2lK - 1) zZ- 1 (72) 

which supports a nonzero magnetization directed along 
the exchange anisotropy axis: 

i '<- 

(73) 

The (+) channel remains gapless. Together with the fi- 
nite spontaneous magnetization this circumstance makes 
the longitudinal spin correlations algebraic and incom- 
mensurate: 

(S?(*)S?, 2 (0)) = (S z ) 2 -f~ 

r_]\x/ ao 
+ const , .„ , cosgoff, (74) 

where qo = S+y/K+/2u. The transverse spin correla- 
tions are short-ranged. 

The ground state is doubly degenerate: the two vacua 
transforming to each other under time reversal. The cor- 
responding topological kinks have a finite mass gap and 
carry the relative spin 



1 



±- 



1 




CSN 



MSN 



D 




FIG. 3: The first path shows one possible way of going 
through CSN, MSN and D phases when z± < z\, by increas- 
ing A and keeping A' to be a constant. The right figure 
shows the qualitative change of the spin current (solid line) 
and the demerization order parameter (dashed line) following 
this path. In the second path by decreasing A' and keeping 
A at constant value, the ladder enters the F phase. 



VII. RG APPROACH AT A-D BOUNDARY 

The results of sections IV and VI are valid far enough 
from the boundary between sectors A and D where one of 
the two twist operators, 0i, 2 , is strongly relevant while 
the other is strongly irrelevant. On the boundary K- = 
1 separating these sectors both twist operators become 
marginal. Therefore we can expect that in the immediate 
vicinity of the boundary, 



K- =1-6. 



(75) 



far away from the SU(2)-symmetric point, i.e. K+ > 1 
in the sense that K + — 1 = 0(1), the infrared behav- 
ior of model will be controlled by the interplay between 
the two parity-breaking operators with a nonzero confor- 
mal spin, 0i ;2 , and the longitudinal (conformal-scalar) 
terms d x ^^d x ^^ responsible for renormalization of the 
coupling constants. So the starting low-energy model 
should therefore contain both twist terms: 



/(+) 



H = Hq 

+ 7i9 K 6 + sin/36- 



7 2 ct.<I> + sin/3$_ 



(76) 



Here 71,2 differ from Ai j2 by some multiplicative factors, 
and 



0- 



47T 
J 



^kKZ 



1 



5_ 



0{6l) 



(77) 



Notice that the two twist terms in (|76|l contain vertex 
operators (the sines) of mutually dual and nonlocal helds, 
0_ and Models of this kind cannot be treated by 
the variational method used in the preceeding sections. 
This is why in this section we address the RG flow of 
this model which will be studied using a mapping of the 
bosonic Hamiltonian l|76|l onto a theory of four interact- 
ing real (Majorana) fermions (c.f. Ref llOri . 

Let us first make all perturbations in l|76|) strictly 
marginal. This can be done by the following rescaling 
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of the fields in the (-) and (+) sectors: 



IK 



K 



9_ -> (l/^/KI) 9_ (78) 

e+ -» (l/v^r) e+ (79) 



The meaning of the first rescaling is transparent: we en- 
force the twist operators in (|76|l to have the scaling di- 
mension 2. This rescaling generates a current-current 
term in the (-) channel. On the other hand, fusing 
the two twist operators, one generates a similar term in 
the (+) sector; that will renormalize the parameter K + , 
K+ — > K+ = K + q + . Below we will set q + = 1— <5 + assum- 
ing that |<5 + | <C 1. The Hamiltonian (|76|) then acquires 
the form: 



[(d a 



(0,9.)' 



2u V 8 s d x $ sR d x $ sL 



s=± 



(80) 



+ 7i 92,6+ sin V47T0_ + 7 2 9 a; < i> + sin V47r$_ 
where 8± satisfy the initial conditions 

<si 0) =6- (81) 



4 0) = 0, 



It is understood that extra factors appearing due to 
rescaling of K + are absorbed into a redefinition of the 
coupling constants 71 and 72. 

The structure of (|80|) immediately suggests mapping 
onto four real (Majorana) fermions, 6 (a = 0,1,2,3). 
This can be done using a correspondence 



(*+,e H 



(6,6), ($_,e_) 



and standard fermionization rules for the currents and 
vertex operators. The resulting theory is given by the 
Euclidean action describing four degenerate massless 
fermions with a chirally asymmetric interaction: 



3 

E 

a=0 



d 2 z {£ a d£ a + £ a d£ a ) 



+2irv J d 2 z [8+ 6666 + <5- 6666 

+ 7+ (6666 + 6666) 
+ 7-(6666 + 6666)]- 



(82) 



Here 6(- 2 an( A 6(^) are holomorphic (left) and anti- 
holomorphic (right) components of the Majorana fields, 
z = vt + \x and z = vt — ix are complex coordinates, 
d = d/dz, d = d/dz, and 



r 3/2 



7± 



o 



2iru 



(71 ±72) 



(83) 



Due to its chiral asymmetry, the interaction in l|82() 
gives rise to renormalization of the velocities already on 
the one-loop level. For this reason we will discriminate 
between the velocities of different Majorana species, and 
set v a = v(l — An p a ) (a = 0, 1, 2, 3), where the dimension- 
less parameters p a are subject to renormalization with 



initial conditions p a 



(0) 



0. 



Using the standard fusion rules for fermion fields^ 3 -, one 
can easily derive the following one- loop RG equations: 



= -27+7-, 
7+ = <*-7-> 

Pi = f>2 
P3 = 7+) 



8- 

7- 



0, 

: 5-7+! 



7+ +7-- 



Po = 7-, 



(84) 
(85) 
(86) 
(87) 



where g = dg(l)/dl, I = \n(L/a). 

First of all, we observe that the coupling constant <5_ 
stays unrenormalized: 

<5_(0=<y-(0) = *_ (88) 
Representing A± as 

7± = 5i ± 52, 5i,2 = (7r 3/2 a/27ra)7i i2 
we rewrite the first, third and fourth RG equations as 



5+ = -2{gf-gl) 
91 = S-gi, gi = -S-g 2 



(89) 
(90) 



We see that, depending on the sign of 8-, either gi(l) or 
<?2(0 grow up upon renormalization: 

(a) 8- > 

gi (l)=g[ 0) e 5 - 1 , g 2 (l)=g^e- s - l ^0 (91) 

Strong-coupling behavior of gi(l) in (|91|l is associated 
with a dynamical generation of a mass gap 



mi cx | 5 i| 1/5 -- 



(92) 



(b) S- < 



m (°) I*- 12 en C°) -15-12 
52(0 =5 2 e 1 1 , 9i(l) =9i e 

Here the mass gap is estimated as 

m 2 oc l^l 1 /!*-!. 



(93) 



(94) 



The cases (a) and (b) describe the CSN and F phases, 
respectively. Estimations H92|) . (|94|l are consistent with 
the power-law scaling of the corresponding mass gaps, 
Eas. H4(jfl and (|72|l . In both cases, flows to strong 

coupling. It goes to large negative values in the case (a), 
implying that K+ becomes even larger upon renormal- 
ization. In the case (b) it flows to large positive values; 
so the effective K + significantly reduces, and that might 
indicate the importance of the neglected twist operator 
O3. Stability of the F phase is therefore under question. 

Exactly at the boundary between sectors A and D 
8- = 0. In this case both g\ and g\ stay unrenormal- 
ized. Moreover, 



8 + (l) = -2{gl-gl)l. 



(95) 
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So, if in addition we set g\ =52, the 1-loop RG will 
display a weak-coupling regime for all coupling con- 
stants. This is the self-dual point of the model where 
the interaction is not renormalized: for all effective cou- 
plings parametrizing interaction the /3-function vanishes. 
Amaizingly, in this case the Majorana action Ij82(l de- 
couples into two chirally asymmetric, independent parts, 
S = 51 + 5*11, where 

Si = J d 2 z (j2 e&c + + 9 eeeA m 

and Su is obtained from Si by reversing the chiralities 
of all the fields. Notice that even though the present 
case corresponds to an essentially anisotropic regime (re- 
member that we are far away from the SU(2)-symmetric 
point of the model), the effective theory on the bound- 
ary between sectors A and D (if_ = 1) with the self- 
duality condition gi = g 2 exhibits an enlarged, chiral 
SO (3) <E> SO (3) symmetry. Consistent with this symme- 
try is renormalization of the velocities. The velocity of 
the singlet fermion, £°, stays intact: po — 0. However, 
the triplet velocity is renormalized. The RG equation 
pt — 4g 2 (pi = p t i = 1, 2, 3) shows that A-np t (l) increases 
upon renormalization and reaches values of the order of 
1 in the region where g 2 l ~ 1. This sets up an infrared 
energy scale in the problem, wo ~ Aexp(— const/g 2 ), 
at which the triplet collective excitations soften signifi- 
cantly. This is in agreement with the exact results for the 
spectrum of model (|96[) . recently obtained by Tsvelik^. 

Interestingly enough, the exact solution^ shows that 
the chiral SO (3) symmetry of the action (|96[) is spon- 
taneously broken at T=0, and the ground state of the 
model represents a "chiral ferromagnet" characterized by 
a nonzero expectation value of the vector current: 

(I) ^ 0, I a = -(i/2)e ab %t; c . 

Similarly, for action 5*11 

(I) ^0, I a = -(i/2)e ab %£ c . 

As long as the actions 5*1 and Su are decoupled, there is 
no correlation between (I) and (I), or equivalently, be- 
tween the magnetization m = (I) + (I) and spin cur- 
rent j = (I) — (I) . Such correlation appears upon devia- 
tion from the A-D boundary since in this case chirally- 
symmetric terms that couple the actions Si and Su (and 
also introduce a finite XXZ anisotropy) are generated. 
Thus, in the A-vicinity of the A-D boundary m z 7^ 0, 
j z = 0, whereas in the D-vicinity the situation is just 
inverted: m z — 0, j z ^ 0. So, the resulting picture 
at the A-D boundary depends on the side from which 
this boundary is approached, implying that the CSN - F 
transition is first-order. 

Even though the action S = Si + 5*11 provides the sim- 
plest field-theoretical model for a frustrated ladder with 
a chirally asymmetric interaction and, hence, is quite in- 
teresting in its own right, we will refrain from its further 



discussion because it does not account for the low-energy 
properties of the zigzag spin- 1/2 ladder with the generic 
SU(2) symmetry (the point K + = K_ = 1). 

VIII. CONCLUSIONS 

In this paper we have analyzed the phase diagram of 
the spin-1/2 anisotropic zigzag ladder with a weak inter- 
chain coupling (J' <C J). Using the Abelian bosonization 
method combined with a variational approach, we have 
found that, depending on the anisotropy parameters, the 
system occurs either in the parity and time-reversal sym- 
metric, spontaneously dimerized phase, or in one of those 
phases in which either parity is spontaneously broken 
while time reversal preserved, or vice versa. These are 
the critical and massive spin nematic phases (CSN, MSN) 
and the critical ferromagnetic (F) phase. We have shown 
that the CSN phase extends well beyond the XX limit 
and covers broad regions A and B in the parameter space 
of the XXZ model (see Figs. 2). Each of these phases 
is characterized by topological excitations carrying frac- 
tional quantum numbers. 

Starting from a vicinity of the XX point, we ad- 
dressed the nature of the transition between the CSN and 
D phases taking place upon increasing the intra-chain 
anisotropy parameter A, say, at a fixed positive value 
of A'. Typical curves are shown in Fig. 2. In these two 
Figures we show two possible scenarios whose realization 
depends on the ratio (zx/z±) between the amplitudes of 
the main competing perturbations - the twist operator 
0i and the dimerization field T>. The reason we consid- 
ered each of these scenarios on equal footing is due to the 
fact that the relationship between the parameters of our 
bosonized model, Eq. l|25|l . and those of the microscopic 
Hamiltonian @ is nonuniversal and, hence, known only 
by the order of magnitude. If z± > z\, the CSN - D 
transition is first order. In the opposite case, z± < z%, 
the CSN and D phases are sandwiched by the MSN phase 
characterized by the coexistence of a finite spin current 
with dimerization. Then the variational approach uan- 
mbiguously shows that CSN-MSN and MSN - D tran- 
sitions are continuous, even though it is inadequate to 
identify their universality classes. We believe that the fi- 
nal choice between the two possibilities discussed in this 
paper will be made in future numerical work (The accu- 
racy of the recent DMRG calculations- 11 ' 12 for the S=l/2 
zigzag ladders was reported to be inadequte to resolve 
this issue. Moreover, the non-linear a model approach of 
Ref. UtI which excludes the massive chiral phase for half 
integer spins, is only valid in the vicinity of the classical 
Liftshitz point j = 1/4.) 

Starting from the region A occupied by the CSN phase, 
one can also keep the in-chain anisotropy intact and vary 
continuously the interchain anisotropy. In particular, one 
can smoothly go from the case of an antiferromagnetic 
interchain coupling [gu > 0) to the case of a ferromag- 
netic coupling (<jni < 0) (see Fig.l). We have shown that 
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in such situation the ladder crosses over from the CSN 
phase to the F phase, the latter being dual to the former. 
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Its minimization yields the following set of equations: 

M K ~ (zj_M+ + sin 7 cos £ =p z iJ+ cos 7) = 0, 
M+ + {z ± M K - cos 7 sin C T £ 3 Q_ cos C) = 0, 
J+TziM K - s\n-i = 0, 
Q-TZ3M++ sinC = 0, 
M_(l T Zi J+K-M*-~ 2 sin.^ 



- zj_if_M+ + M"- A cos 7 cos C) = 0, 
M+(l t ZsQ-K+M? + ~ 2 sinC 



ziK+Mf- M+ + " cos 7 cos C) = (A3) 



rK 



APPENDIX A: MORE ABOUT SECTOR C 

In this Appendix we address the role of the so far ne- 
glected twist perturbation A3O3 which becomes relevant 
in sector C. Adding this term to the effective Hamiltonian 
leads us to the following theory: 



(+) 




(-) 




K 

— cos J^kK~<^ + cos J^kKZQ- 
2Ai 



-d x Q + sin 



-.d x §- sin y/4:irK + $ + . 



(Al) 



It is convenient to introduce dimensionless notations for 
the coupling constant 



Z3 



TV X^a 

K- v 



and the z-component of the relative spin density, 

The variational approach we followed in section V is 
straightforwardly generalized for the present case. As 
compared to section V, here we have two additional vari- 
ational parameters: the relative spin density Q_ and a 
mixing angle Q for the field <!>+. The variational energy 
then depends on six variables: 



£c = \{M 2 + 



M 2 _+Jl 



Q 2 _) 
IM 

cos 7 cos £. 



=F zi J+ M K sin 7 =p z 3 QM+ + sin ( 
- z ± M++M*- 



There exist solutions of these equations in which the 
second twist perturbation A3O3 plays no role: 

(i) 7 = C = , J+ = Q- = - D phase; 

(ii) 7 = vr/2 , C = , Q- = - CSN phase;; 

(iii) 0<7<tt/2,C = 0,Q_=0 - MSN phase. 

There exists a pair of solutions which are "dual" to (ii) 
and (iii), i.e. can be obtained from the latter by the 
replacements z\ — > Z3, J + — > Q_, M + «-> M_: 

(iv) 7 = , ( = ir/2 , J + = 0. This is a critical phase 
with a nonzero relative magnetization (CRM), Q_ ^ 0; 

(v) 7 = , < C < tt/ 2 , J+ = - the massive 
version of the above phase (MRM). In these two phases 
the twist operator 0\ plays no role. 

There also exist solutions in which both twist pertur- 
bations are effective. One of them corresponds to the 
case 

(vi) 7 = 7r/2 , £ = 7r/2 , in which the z^-perturbation 
is ineffective and the variational energy decouples into 
a direct sum £csn + &CRM- The resulting phase is 
fully gapped and represents a mixture of CSN and CRM 
phases - mixed (M) phase with nonzero J+ and Q_. 

The case of arbitrary values of the mixing angles, 
7, £ 7^ 0, 7r/2, should be abandoned because, as follows 
from Eas. (|A3|l . it requires that z\z\ = z\, - a condition 
which represents just a point in the parameter space of 
the model and which, on the other hand, cannot be sat- 
isfied for all coupling constants being of the same order. 

The minimal value of the variational energy is again 
given by Eq. Il-Sl'l) . From this expression it is obvious that 
in sector C (K± < 1) the M phase has a lower energy than 
each of its "constituents" , i.e. CSN and CRM phases. So 
we are left to find out if the M and MRM phases can 
compete with the D and MSN phase. 

Consider the MRM phase assuming the most favorable 
condition z% > z±. Since the MRM phase is "dual" to 
the MSN phase, from Ea. (|68|l we can read off the range 
where it can exist: 



1 - (1 - K+) 



ln(z 3 /z ± 
ln(l/z 3 ) 



< < 1. 



(A4) 



We see that, except for an extremely unrealistic case 
z\ > z±, the condition 
(A2) negative semiaxis A"_ = 



4} determines a vicinity of the 
1, K + < 1, which is located in 
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the unphysical part of the (K + , K_) plane, well beyond 
sector C. Thus the MRM phase should be abandoned. 

Let us compare the energies of the M and MSN phases. 
In both cases the mass M_ is given by the same ex- 
pression, so we only need to compare the masses in the 
(+) channel. Comparing the mass M+ in the M phase, 

M+ ~ z 3 1/(1 ~ K+) , with that in the MSN phase, Eq. j52). 
we find that, except for extremely small values of Zj_, 
namely zj_ < Z\Z%, the MSN phase is always more favor- 
able. 

Finally, we are left to compare the energies of the M 



and D phases. On one hand, in the D phase we are 

i i 

below the line l|t)8l) . This means that zl_ D > z 1 ~ 
implying that the mass gap of the D phase is greater 
than the mass M_ of the M phase. On the other hand, 

to the left of the MSD-D transition line (|A4|) we have 

i i 

the condition z^ D > z 3 + that tells us that the mass 
of the D phase is greater than the mass M+ of the M 
phase. Consequently, the D phase is energetically more 
favorable than the M phase. 
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